Segregation of granular binary mixtures by a ratchet mechanism 
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We report on a segregation scheme for granular binary mixtures, where the segregation is performed 
by a ratchet mechanism realized by a vertically shaken asymmetric sawtooth-shaped base in a quasi- 
two-dimensional box. We have studied this system by computer simulations and found that most 
binary mixtures can be segregated using an appropriately chosen ratchet, even when the particles 
in the two components have the same size, and differ only in their normal restitution coefficient or 
friction coefficient. These results suggest that the components of otherwise non-segregating granular 
mixtures may be separated using our method. 

PACS numbers: 45.70.M, 05.60.C, 05.40. J 

While segregation is often an undesired effect, some- 
times separating the components of a granular mixture 
is the ultimate goal. Since ancient times sieves have been 
used by humans to separate small grains from bigger 
ones. But nature manages to separate different kinds of 
grains also without sieves. Many of the segregation pro- 
cesses in granular matter [0-^ have recently been studied 
in great detail such as segregation according to particle 
size, shape or friction properties in a shaken box ^-|^ , in 
a rotating drum [Q- or when poured into a thin box 
pTj-p^. However, so far no method has been proposed 
for segregating particles which only differ in hardness. 

The setup we use for segregation is as follows: particles 
of a binary mixture are falling into a two-dimensional box 
from above at a specified place. The base of the box hav- 
ing an asymmetric sawtooth profile is shaken harmoni- 
cally in vertical direction. The transport properties of ho- 
mogeneous granular media in a similar setup were studied 
earlier p^ ] inspired by recent progress in the theoretical 
understanding of molecular motors . In the corre- 

sponding models, known as thermal ratchets, fiuctuation- 
driven transport phenomena can be interpreted in terms 
of overdamped Brownian particles moving through a pe- 
riodic but asymmetric (typically sawtooth shaped) po- 
tential in the presence of nonequilibrium fluctuating 
forces (such as periodic driving or switching between po- 
tentials). When a particle jumps out of the box at either 
the right or the left boundary, it is removed and counted, 
and finally the segregation quality is determined. For de- 
tails on the geometry see Fig. ^ The segregation quality 
isq= (2max{7Vi<_-h7V2^,iVi^-hiV2^}/7V-l) x 100%, 
where i^i^) denotes the number of particles of 

the zth component {i ~ 1,2) leaving the box on the left 
(right) side, and N is the number of all particles in the 
mixture. Thus the quality of random segregation, when 
the particles go to the left or right side with equal prob- 
abilities, is 0%. In this Letter we show results for binary 
mixtures in which the components contain equal num- 
ber of particles, in other cases different quality defini- 
tions may also be appropriate. Unlike other segregation 



FIG. 1. Schematic drawing of the segregating setup. The 
particles are falling into the box of width L at aL from the 
left side of the box (0 < a < 1) from height H with load rate 
R. The base is shaken sinusoidally with amplitude A and fre- 
quency /, and may be tilted by angle 9, which is positive if 
the left end is higher than the right one. The particles within 
a mixture component are exactly the same if not stated other- 
wise. Inset: The shape of one sawtooth is described by three 
parameters: width w, height h and asymmetry parameter a. 



phenomena, in which segregation is due to the collective 
behaviour of the grains, here the interaction between the 
base and the individual particles is dominant, and the 
efficiency of segregation usually decreases with increas- 
ing number of particle-particle collisions. However, set- 
ting a sufficiently low load rate R, the quality can be 
high, and still, "parallelizing" the procedure the separa- 
tion capacity can be large (see bellow) . In the simulation 
an event-driven algorithm is applied with a hard- 
sphere collision model |l^, in which the sphere-shaped 
particles have five parameters: mass m, radius r, normal 
restitution coefficient e, friction coefficient /i and maxi- 
mum tangential restitution coefficient f3o. The particles 
can rotate around the axis going through their center and 
perpendicular to the plane of the box, their moment of 
inertia is 2/5 mr^ about their center. Since the mass of 
the particles does not play a role in collisions with the 
base but only in binary collisions, it is enough to spec- 
ify that the particles have the same mass density, so the 
mass is cubically proportional to the radius. 
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If there are only few particles in the box at the same 
time, then, as a first approximation, the interaction be- 
tween the particles can be neglected. Therefore, we in- 
vestigated the motion of one particle in an infinitely wide 
box in detail, and found it to be chaotic in most 
cases. A similar, but simpler model was also reported to 
show chaotic behaviour [ pO[ . Depending on the param- 
eters, it is possible that the particle follows a periodic 
trajectory, travelling with velocity v — fwb/c, where b 
and c are integer numbers, meaning that in one period, 
which lasts for c vibration cycles, the particle jumps over 
b teeth. These periodic trajectories are not interesting 
for practical applications for the following reasons: (1) 
the transients are usually very long, and therefore may 
be more important than the asymptotic periodic trajec- 
tory for the segregation behaviour, (2) periodic trajecto- 
ries are not robust against collisions with other particles 
and other sources of noise, and, last but not least, (3) 
for certain conditions two periodic trajectories with op- 
posite directions can coexist for the same type of parti- 
cles. In this case one cannot predict on which side the 
particle tends to leave the box. We explain below how 
these periodic trajectories can be avoided when search- 
ing for the ratchet parameters suitable for segregating a 
given binary mixture. In the chaotic regime, however, 
the time evolution of the particle's horizontal position 
can be well described as drift-diffusion. The connection 
between chaotic motion and diffusion has been investi- 
gated extensively recently pl|-|23[| . A simple explanation 
for this drift-diffusion behavior can be that on time scales 
larger than the typical time it takes for the particle to 
jump to another sawtooth, the kicks of the base can be 
considered to be independent of each other. Further- 
more, the asymmetry of the ratchet leads to an average 
velocity in the left or right direction. Consequently, the 
horizontal motion in the chaotic regime can be described 
statistically by two parameters: a drift velocity v and a 
diffusion coefficient D. An example for the drift-diffusion 
motion can be seen in Fig. 

The observation that the horizontal motion can be 
well approximated by drift-diffusion enables us to predict 
segregation-related quantities . First of all, the prob- 
ability that a particle jumps out of the box through the 
right (n^) or the left {n^) boundary is: n^(u, L, a) = 
^^^^-Lu and n^(u, L, a) = 1 — n^{u, L, a), where we 
introduce the notation u — v/D, since the probabili- 
ties depend on the drift velocity and the diffusion co- 
efficient only through this combination. The asymptotic 
behaviour of these probabilities in L is exponential with 



characteristic length a 



which means that if the 



drift velocity is, say, positive and the box width is large 
enough, most or all of the particles arrive at the right 
end. The approximate explanation for this is the fol- 
lowing: the displacement of the average particle position 
increases linearly in time, while the width of the probabil- 
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FIG. 2. Drift-diffusion of a single particle in an infinitely 

wide system (realized by periodic boundary condition). The 
probability distribution of the horizontal displacement x in 
time t — 2 s (square), t — 10 s (circle) and t = 20 s (tri- 
angle). The dashed curves show the corresponding Gaussian 
distributions with center vt and dispersion \/2Dt as theoreti- 
cal prediction |2^] . The drift velocity v and diffusion constant 
D are determined by line fitting: {x{t)) = vt and o-{t)^ = 2Dt 
(see inset), where <T(t)^ — {x{t)'^) — {x{t))^ . The fitted val- 



ues are v = 1.38 cms and D — 5.38cm s 



The particle 



parameters are r — 1.5 mm, e = 0.6, — 0.3 and /3o = 0.4, 
the ratchet parameters are w = 12 mm, h — 8 mm, a — 0.2, 
e ^0, A = 2mm and / = 20s"^ 



ity distribution is proportional only to the square root of 
time (although this is exactly true only if the box width is 
infinite, for large box widths it is still a good approxima- 
tion). Therefore, on length scales larger than the 
drift dominates over diffusion, so that most likely the par- 
ticle leaves the system at the end towards which it drifted. 
As a consequence, for the segregation of a binary mixture 
the theory suggests that if the drift velocities of the parti- 
cles of the two components (vi and V2) have opposite di- 
rections, then one can obtain an arbitrarily good segrega- 
tion quality by choosing the box wide enough. For a fixed 
box width, the best segregation quality is given by gopt = 
[n^{ui,L,aopt) + n^{u2,L,aopt) - 1] x 100%, where 

Ml < U < U2, and a^pt - t (u2-ui)L 

gives the optimal place of loading to obtain the best pos- 
sible quality. The asymptotic behaviour of gopt as a func- 
tion of L shows that it exponentially approaches 100% 
with characteristic length -I- \u2\) |wiU2|~"'^. Figure 
^ shows an example for segregation of particles differing 
only in friction coefficient. In this example the parti- 
cle parameters are: ri = r2 = 1.5 mm, ei = 62 = 0.4, 
/ii = 0.1, /i2 = 0.3, /3oi = /5o2 = 0.4, the ratchet param- 
eters are: w — 8 mm, h — 6 mm, a ~ 0.4, 9 — —0.08, 
A = 2mm, / — 20 s^^, and the segregation parameters 
are: a = 0.51, H = 1cm, R — 0.5 s~^. The correspond- 
ing diffusion parameters are vi — — 1.58cms"^, Di = 



1.36 cm^ s 



V2 



1.86cms"i and Do = 1.83 cm^s^^. 
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FIG. 3. Segregation of particles differing only in friction 
coefficient. The segregation quality (solid line) rapidly grows 
to 100% with increasing box width, the dashed line shows 
the theoretical prediction. At small box widths the devia- 
tion from the theoretical prediction (dashed line) is due to 
that (1) the drift-diffusion approximation is valid on length 
scales larger than the sawtooth width, and (2) the starting 
state (position and velocity) of a particle is untypical to the 
drift-diffusion motion. Inset: The drift velocity and diffusion 
coefficient (measured in an infinitely wide box) as functions 
of the friction coefficient /i. The drift velocity changes its sign 
twice, at ^ = 0.14 and at /i = 0.4, and is only coincidentally 
zero at /i = 0. 



the characteristic length is 1.85 cm. The results are ob- 
tained from the segregation of 20,000 particles, 10,000 in 
both components. 

The ratchet segregating a certain binary mixture most 
efficiently is searched in the following way: for both com- 
ponents the diffusion parameters of the one particle mo- 
tion (wi, Di and V2, D2) are measured using many differ- 
ent ratchets. Then we select those ratchets for which the 
drift velocities have opposite directions (it may happen 
that no such ratchet is found). Setting the box width 
to a reasonable value, for each of the selected ratchets 
the best segregation quality is predicted, and the ratchet 
with the highest segregation quality is chosen for segre- 
gating the mixture. However, it is possible that with this 
ratchet one or both of the particles have periodic trajec- 
tories, which is undesired for segregation. We describe 
here one possible solution to this problem: any kind of 
noise can destroy periodic trajectories. For example, in 
the computer simulation the angle of the relative velocity 
after a particle-base collision is changed by an amount 
uniformly chosen from an interval [—(50,50]. We found 
that a rebounding angle noise Sip ~ 0.05 (measured in ra- 
dian) is enough to destroy the periodic trajectories, dras- 
tically changing the drift velocity. Hence, the diffusion 
parameters for both particles and each ratchet are also 
measured with 54) — 0.05 and 0.1. Then, for selecting 
a ratchet it is not enough that the drift velocities have 
opposite direction when 5(f) = 0, but their values should 
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TABLE I. The parameter values used to produce 3'^ = 27 
particle types and 5^ = 3125 ratchet types. 



not change too much when the noise level is increased to 
0.05 and 0.1 (we allowed a maximum relative change of 
50% and 70%, respectively). The best ratchet is chosen 
from this restricted set by a rule which takes into account 
that (1) the quality should be high also in the noisy case, 
(2) \vi\ and |w2| should be as large as possible to allow a 
high load rate, (3) aopt should not change much in the 
presence of noise, otherwise the segregation quality may 
decrease much due to interaction between particles. This 
choice of noise is not only good for avoiding periodic tra- 
jectories, but also may serve as a first approximation for 
taking into account the deviation of the particles' shape 
from a sphere. Therefore, the ratchets selected by this 
method are robust to some extent against deviations in 
shape and presumably in other parameters as well. 

We checked the practicability of this procedure by find- 
ing appropriate ratchets to segregate {^^) = 351 different 
binary mixtures composed from 27 particle types, mea- 
suring the diffusion parameters for each type by simulat- 
ing 3000 seconds of particle trajectories for 3125 ratchets. 
The parameters which were varied are listed in Table |, 
the maximum tangential restitution coefficient and the 
vibration amplitude were fixed /3o — 0.4 and A = 2 mm, 
respectively. Then we performed a segregation simula- 
tion for each mixture with the selected ratchet (for three 
mixtures no proper ratchet was found in this set). The 
results show that in most cases good segregation quality 
can be achieved even if the rebounding angle noise is rel- 
atively high or the particles within the components are 
not uniform (see Fig. |^). 

If no ratchet is found, even in an extended set, which 
segregates a certain mixture (i.e., the drift velocities al- 
ways have the same sign), then there is an alternative 
procedure which we describe here only briefly. Unless the 
particles in the components are exactly the same, there 
exists a ratchet for which ui and U2 are different. If, 
6.17., < Ml < U2^ then choosing an appropriate loading 
place near the left end, more particles from component 1 
will be collected on the left side than from component 2. 
Then the particles collected on the right side are reloaded 
again and again. For a given box width and number of 
reloads one can determine the optimal load place to ob- 
tain the best possible segregation quality, which tends 
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Segregation quality, q [%] 

FIG. 4. Segregation quality for 351 binary mixtures. In all 
cases the segregation parameters are L = 30 cm, H = 3 cm 
and R = 0.25 s~^. a, The particles are uniform within the 
components, and the rebounding angle noise is Sep — 0. The 
number of mixtures in the last column (quality is better than 
95%) is 219, and the quality is better than 70% for 320 mix- 
tures, b, The same as in a but the rebounding angle noise is 
5(j) = 0.3. The quality is better than 70% for 270 mixtures, c, 
The same as in a but the particles are not uniform within the 
components, all of the parameters (r, e, fi and /3o) are varied 
by maximum 10%. The quality is better than 70% for 263 
mixtures, d, The rebounding angle noise is Sep = 0.3 and the 
particle parameters are varied by maximum 10%. The quality 
is better than 70% for 220 mixtures. 



to 100% with increasing number of reloads (results not 
shown here). 

One may think that the method presented here is not 
efficient for segregating real granular mixtures, since only 
a few grains can be present in the box at the same time, 
allowing only a low load rate. However, it is very easy to 
"parallelize" the procedure: as the box is essentially two- 
dimensional, many boxes can be placed onto a shaking 
machine, and possibly the walls between the boxes can 
be omitted. A rough estimation shows that the capacity 
of such a machine would be comparable to that of the 
machines in use today in the industry for cleaning, e.g., 
cereal grains (in which case the capacaity is in the order 
of 1 ton per hour). We are planning to carry out ex- 
periments to check if mixtures containing particles of the 
same size differing only in normal restitiution coefficient 
or friction coefficient can be segregated by this method. 

We presented a computer simulation study of a method 
for segregating a binary granular mixture. The segrega- 
tion is performed by a ratchet mechanism, and in contrast 
to other segregation schemes, here not the collective be- 
havior of the particles is dominant but the interaction 
between the base and the individual particles. We found 
that good segregation quality can be achieved even if the 
particles of the two components differ only in friction co- 
efficient or hardness. 
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